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Abstract 

Any quantum system with a non-trivial Hamiltonian is able to simulate any other 
Hamiltonian evolution provided that a sufficiently large group of unitary control op- 
erations is available. We show that there exist finite groups with this property and 
present a sufficient condition in terms of group characters. We give examples of such 
groups in dimension 2 and 3. Furthermore, we show that it is possible to simulate an 
arbitrary bipartite interaction by a given one using such groups acting locally on the 
subsystems. 

1 Introduction 

Simulating Hamiltonian evolutions of arbitrary quantum systems is an interesting ap- 
plication for future quantum computers. Historically, the idea of simulating Hamilto- 
nian time evolutions was the first motivation for quantum computation Q. Whereas 
in early works on this problem the desired Hamiltonian was proposed to be simulated 



by a discrete sequence of gate operations (see e. g. [23]), a more control theoretic for- 
mulation of the problem has become popular recently II' 1> H> 0- In this 
formulation one assumes that the dynamics of the quantum computer is determined 
by its Hamiltonian together with external control possibilities. Here the task is to 
simulate an evolution that would have occurred under some other Hamiltonian by in- 
terspersing the natural time evolution with control operations. More explicitly, one 
assumes that the natural evolution exp(-iHt) alternates with fast implementations of 
unitary operations Vj, i.e., the resulting evolution is given by 

exp(—iff rjvi)V/v ■ ■ ■ exp(—iHT2t)V2 exjp(—iHTit)Vi , (1) 
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where the relative times between the control operations are given by r, for i = 1, . . . , N. 
If the t is small compared to the time scale of the evolution exp(-iHt), the resulting 
evolution is approximatively given by 

exp(-iHt) . 

Here H is the average Hamiltonian 

N 



i=l 

where we set U{ := YYj=iVj- Important characteristics of simulation schemes like (||) 
are the amount of operations performed (complexity) and the total operation time 
(overhead). The average Hamiltonian method is the basis for simulating Hamiltonians 
by a given one and has applications which go beyond the usual aims of quantum 
computation. As examples we mention the following applications: 

1. Decoupling / Suppression of Decoherence 

The time evolution on the joint Hilbert space of a system and its environment 
can be interspersed with transformations on the system's space alone in such a 
way that the net effect is a separate ( "decoupled" ) time evolution of the system 



and the bath [35]. This is a generalization of decoupling techniques in Nuclear 
Magnetic Resonance B. 

2. Time inversion 

The natural time evolution t t— > exp(-iHt) with t > can be conjugated by 
unitary transformations in such a way that the total effect is a transformation of 
the form exp(-iHs) with negative s. Remarkably, there are schemes for inverting 
unknown Hamiltonians. This fact is closely related to the existence of decoupling 



schemes [35]. Time inversion for unknown Hamiltonians is a useful primitive in 



quantum process tomography, i.e., procedures that distinguish between unknown 
time evolutions [ 14 ] . The essential idea is that time inversion enables to implement 



transformations of the form exp([H, A]) for arbitrary self-adjoint A even if H is 
unknown. 

3. Generating arbitrary time evolutions with a finite control group 

Assume that the only external control operations of a quantum system are im- 
plementations of unitary operations taken from a finite group S. If the natural 
time evolution is non-trivial and 5 in a sense is large enough, than every unitary 
transformation can be generated by concatenations of the natural evolution and 
elements of S. In particular, this is possible if the natural Hamiltonian H can be 
transformed into any other in the sense of the average Hamiltonian method. 
This is a special instance of the following problem: Given an arbitrary M-linear 
map L on the set of self-adjoint traceless operators, the task is to find a sequence 
of unitary operations such that the system evolves according to the Hamiltonian 
L{H) if its natural (unknown) Hamiltonian H is present. The problems of invert- 
ing or switching off unknown Hamiltonians are special cases of this problem with 
L = or L(H) = —H. In its full generality, procedures for simulating L{H) if an 
unknown H is present can be used as primitives in quantum process tomography 
[14] and quantum control procedures |fl6| |. 
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4. Simulating interactions by other ones 

Consider a bipartite system with Hilbert space C d <S> C d . Assume that the in- 
teraction between both parts is a fixed Hamiltonian on the joint Hilbert space 
which cannot be controlled at all. The only possibilities to control the system 
are implementations of local transformations on each of the subsystems. These 
transformations can be used to imitate other interactions j|, 31]. 
In n-partite systems (C rf )® n , where the interaction between each pair of subsys- 
tems is assumed to be fixed, one can simulate other pair- interactions in the sense 
of the average Hamiltonian method. Time optimization of these simulations leads 
to interesting problems of parallel execution. For n-qubit networks, this problem 
has been addressed in fl34L FL5L EM E2l . 



The paper is organized as follows. In Section || we define what it means to simulate a 
Hamiltonian and explain the physical meaning in the context of average Hamiltonian 
theory §, 0§. 



In Section || we introduce the concept of an annihilator for ci-dimensional quantum 
systems characterizing control procedures that switch off the possibly unknown dynam- 
ics of the system. These procedures directly give decoupling and inversion schemes. 
Explicitly, they can be constructed using nice error bases [ll| which yield minimal 
annihilator procedures attaining the lower bound d 2 on the complexity. We show a 
lower bound on the overhead of inversion procedures of a general, possibly unknown 
Hamiltonian to be d — 1 and an upper bound of d 2 — 1. Furthermore, a lower bound 
on the complexity is shown to be d 2 — 1. 

In Section |4] we address the question how to simulate an arbitrary Hamiltonian by 
any other if only a restricted set of control operations is available. The condition on the 
set of available control operations for making universal simulation possible is stronger 
than the requirement to make annihilation possible. A sufficient condition for a group 
of available control operations allowing universal simulation can be formulated in terms 
of group characters. Groups satisfying this criterion will be called transformer groups. 
We present transformer groups for dimensions 2 and 3. Furthermore, we show in 
Section || that the transformer groups allow to transform any interaction in a bipartite 
system into any other by operating on the subsystems only. 



2 Dynamical control 



For the implementation of a quantum computer it is necessary to control the time evo- 
lution of the used physical system in a universal way. In many physical systems the only 
directly accessible control possibilities are given by a set of control unitaries and the 
system Hamiltonian that cannot be switched off. By applying the control operations 
we effectively change the Hamiltonian into a piecewise constant time-dependent Hamil- 
tonian. The formalism of average Hamiltonian theory [||, 11, 27 allows to solve for the 
resulting time evolution at a time t by writing the evolution of a time independent 
average Hamiltonian H. Following & O, ^7]] we briefly sketch average Hamiltonian 
theory: the overall dynamic after a period of evolution is given by 

U(t) = TexpH / drH(r)) = exp(-iffi) , 
J o 
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where T denotes the Dyson time ordering operator. A solution of this equation is a 
time independent Hamiltonian that would result in the same unitary if it were applied 
over the same period. If the Hamiltonian H(t) commutes with itself at all times we 
have H = f Q l drH(r). However, this is rarely the case. For sufficiently small t, the 
Magnus expansion provides a formal means of calculating the average Hamiltonian: 

H = + fl-W + jjW + . . . (2) 
where the operators H^\H^-\ . . . are the average Hamiltonians of increasing order 

H<® = - tdTH(T), (3) 
T Jo 

BW = ^J\r' j\r"[H(r'),H(r")]. (4) 

We have < 1 and \\H<ti\\ < t/2 since \\[H(t'), H(t")]\\ < 2||#(t , )||||#(t")|| = 2 

and we are integrating over the simplex of area t 2 /2. The norm of the higher order 
terms is bounded by higher orders of t. Therefore for sufficiently small time t the 
resulting unitary U(t) is essentially determined by H^°\ 

We assume that the only directly accessible control possibilities are the unitaries 
in the control set C C SU (d) and assume in addition that they can be performed 
arbitrarily fast compared to the natural time evolution. This socalled fast control limit 
is justified e. g. in NMR because the coupled and local evolutions act on significantly 
different time scales. We will now subject the system to a cyclic pulse train. The pulses 
are assumed to be infinitely short. A sequence consisting of N pulses will be denoted 
by 

P:=(^i,n,... ,V n ,t n ) (5) 

where Vi G C and n > are relative times, i.e. r, > and J2iLi r « = 1- ^ ne pulses 
are applied from left to right. If we apply the sequence over the time t the resulting 
unitary is given by 

N 

U P (t) = Y[exp(-iHTit)V t . (6) 

i=i 

The Tj specify the fraction of time between the pulses. For a cyclic sequence (defined 
by ni=i Vi = 1) we can express the resulting unitary as 

N 

U P {t) = X\U} exp(-iHTit)Ui , (7) 
i=i 

where U { = Y[ l j=1 Vj. Using the identity Z7+ ex.p(H)U = exp(U*HU) wc get 

N 

Up(t) = Y[exp{-iHiTit), (8) 

i=i 
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where Hi = UjHUi. These operators are the Hamiltonians in the so-called "toggling 
frame". Let Adc(H) denote set of conjugates of H 

Ad c (H) = {Adu(H) = U ] HU \ U e C} . (9) 

Then the unitary Up{t) is the solution of a time-dependent Schrodinger equation with 
piecewise constant Hamiltonians in Adc(H). 

The previous discussion motivates the following definition of the notion of simulat- 
ing a Hamiltonian by another Hamiltonian. 

Definition 1 (First order simulation) Let H be any Hamiltonian. We say H can 
be simulated by H with overhead 1, written H -< H , if and only if there are Tj > 
summing up to 1 and Hi S Adc(H) such that H = J2j T iHi, i. e. H can be written as 
a convex combination of conjugates of H by elements of C. H can be simulated by H 
with overhead r iff H -< tH . 

Using this definition of simulation the problem of time-optimal simulation of a 
Hamiltonian is reduced to a convex optimization problem (this has been noted in 

BHD- 



3 Annihilators 

In this section we introduce the concept of an annihilator for a d-dimensional quantum 
system characterizing control procedures for switching off the possibly unknown dy- 
namics of the system. These procedures directly give decoupling and inversion schemes. 
We prove some optimality properties of annihilators and show how a minimal annihi- 
lator can be explicitly constructed using nice error basis. 

Definition 2 (Annihilator) Let P := (Vj.,Ti, V%,t%, . . . ,Vn,t~n) be a cyclic pulse 
sequence of length N, i. e. T{ > 0, Yli=i r « = 1> an ^ YiiLi V{ = 1. Set Ui = YYj=i Vj- 
The sequence P is called an annihilator of dimension d and length N iff 

N 

^nfj\aUi = Q ( 10 ) 
1=1 

for all a £ su(d). An annihilator is called minimal if there is no shorter annihilator. 



Theorem 1 (Minimal annihilator) A minimal annihilator has length d 2 . Further- 
more all relative times Tj of a minimal annihilator are equal. 

Proof: Let (Vi, tl, V2, T2, • • • , Vzv, tjv) be an arbitrary annihilator. The corresponding 
Ui = II}=i an d T i define a realization of the depolarizing channel A on C d by random 



external fields [36] since 

N 



A(p) = J2riU} P U i = l/d. (11) 



i=l 
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This permits to show that N > d 2 as follows: by sending one part of a maximally 
entangled state in <C d ® C d through the depolarizing channel we end up with 

the maximally mixed state 1/d 2 . Therefore we need at least d 2 unitaries since the 
rank of each (1 (g> uj)\^)(^>\(l &> U{) is 1 and they must sum up to d 2 (the rank of the 
maximally mixed state). For N = d 2 all Tj must be equal for entropy reasons (see [M |, 
page 518) since for r, / jp 

S (jC 7 "^ 1 ® ® ^) j < ^(n,r 2 , . . . ,T0) < log 2 d 2 = S(l/d 2 ) , 

where 5 denotes here the von-Neumann and H the Shannon entropy. One proof of 
the existence of a minimal annihilator is based on the concept of nice error basis to be 
introduced in the next subsection. □ 
Theorem [l] shows that a minimal annihilator corresponds to a unitary error depo- 
larizer [^] and thus to a unitary basis. 

3.1 Nice error bases 

In this section we deal with the problem to construct a minimal annihilator. For this 
we construct bases for the vector-space C dx d of d x d-matrices which consist entirely of 
unitary matrices and are orthogonal with respect to the trace inner product. One way 



of constructing such bases relies on the concept of nice error bases. We refer to [2C] and 



Tql for an overview of this method and mention that nice error bases are used in the 



construction of quantum error control codes [||, [R], IS, 2£]. They are also of interest 
in the theory of noiseless subsystems [21, 35] and in connection with the development 
of quantum authentication codes |ffl| . 



Definition 3 Let G be a group of order d with identity element e. A nice error basis 
on C d is a set £ = {U g E U{d) \ g E G} of unitary matrices such that 

(i) U e is the identity matrix, 

(ii) tr U g = d5 9i e for all g G G, 

(in) U g U h = a(g, h)U gh for all g,h £ G, 

where the factor system a(g, h) is a function from G x G to the set C x := C \ {0}. 



In |18[ it was shown that the map g ^ U g defines a projective representation of G; 
this is a consequence of conditions (i) and (iii) . Condition (ii) shows that the matrices 
U g are pairwise orthogonal with respect to the trace inner product (A, B) := tv(A^B) / d. 
Hence, a nice error basis is an irreducible projective representation of the finite group 
G. The group G itself is also called index group since its group elements index the 
elements of the nice error basis £. 

Note that in general the group generated by the matrices U g for g G G will be larger 
than G, since these matrices are not closed under multiplication. A well-known theorem 
from projective representation theory (cf. Jl^, Theorem V.24.6], [jl3| , Theorem 11.15]) 
states that it is always possible to switch to an equivalent projective representation such 
that the images U g generate a finite group G (see also |l^]). This group is called the 
abstract error group corresponding to £. Whereas g *— > U g is an irreducible projective 
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representation of G, this yields an irreducible ordinary representation of G. It is a 
well-known fact that G is a central extension of G (cf. [|nj): denoting the center of G 
by C(G) this means that G/((G) = G. 

Given a nice error basis {U g \g £ G}, then the abstract error group is isomorphic 
to the group generated by the matrices U g . The assumption that the factor system a 
is of finite order ensures that the abstract error group is finite. 

Example 1 (Heisenberg group) The discrete Fourier transform of length d 6 N is 



the unitary transformation defined by DFT^ := -^=\ 

the primitive <i-th root of unity e 2m / d . Define 
0, . . . , d — 1}, where 



k-l 



s 



/0 1 



V 1 



\ 

1 



)k,l=0,... 

( 1 



_i, where u> denotes 
= 0,... ,d- l,j = 



T : = DFTT • S ■ BFT d 



V 



Here again uj is the primitive d-th root of unity e 27n,//rf . Then Ed is a nice error basis on 
C d showing the existence of nice error bases for any dimension dGNof the underlying 
system. The index group in this case is the abelian group G = x whereas 
the corresponding abstract error group is a nonabelian group isomorphic to a semi- 
direct product G = (Zd x Z£) x Z& (the so-called Heisenberg group). The projective 
representation of G leading to the error basis is defined by mapping the generators 
of G as follows: (1, 0) i— > S and (0, 1) i— > T. The identity ST = ujTS is readily verified 
which shows that the commutator subgroup of G is contained in the center Q{G). This 
also shows that the factor system a corresponding to the projective representation of 
G defined £4 is given by 



<*((i,j), (k,l)) = u 



for all (i,j),{k,l) £ G. 



We give a brief account of some general properties of nice error bases (see also |l~8|1). 
A complete classification of abstract error groups on C d for 1 < d < 11 was given m 
fl8|| . Index groups of abstract error groups are in general not abelian: in 18 ] a family of 



groups having nonabelian index groups was constructed. It is known that all abstract 
error groups are solvable. Moreover, it is known that all solvable groups can occur as 
subgroups of index groups of nice error bases. On the other hand, it is known that not 
all solvable groups can occur as index groups. 



3.2 Averaging and Annihilation 

Using the concept of abstract error groups we describe the idea of switching off an 
interaction by averaging over a group. Whereas usual techniques are based on ordi- 
nary irreducible representations [j35f |, the following lemma shows that averaging over a 
projective irreducible representation also projects onto the set of scalar matrices. 
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Lemma 1 Let M G C , G be a finite group, and R : g >—> U g £ U(d) an irreducible 
projective representation of G. Then the following equation holds: 



1 

W\ 



E u l MU 9 



g&G 



tr(M) 
d 



Proof: We have seen that each projective representation of an index group G with 
associated factor system a gives rise to an ordinary representation of the corresponding 
abstract error group G and that G is a central extension of G. It follows that {U g : 
g G G} is a set of coset representatives for C(G) in G, i.e., 



G 



U C(G)U g . 

geG 



Each element a G G has a unique factorization of the form a = zg where z G C(G) and 
g £ G. From Schur's Lemma (cf. ||, Section 2.2]) follows that for M G C dxd , G a 
finite group, and R : a ^ U a € U(d) an irreducible (ordinary) representation of G the 
following identity holds: 



|G| ^ 



tr(M) 



(12) 



7£G 



Using this we obtain 



I ^ UtMU c 



1 1 creG 

tr(M) 



d 



1, 



where the last line is due to Schur's Lemma (^) for ordinary representations. □ 



3.3 Decoupling 

We consider a bipartite quantum system (e. g. a system coupled to a bath) living on 
the joint Hilbert space Us ®H.b- Here Tis and Hb denote the Hilbert spaces of S 
and B respectively. Let su(Ti) denote the Lie algebra of traceless self-adjoint matrices 
acting on the Hilbert space TC. The Hamiltonian can be written as 

H = H s ®1b + 1s®H b + H S b, (13) 



where H$ G suiTLs) is the free system Hamiltonian, Hb G suijis) is the free bath 
Hamiltonian, and H$b describes the coupling between the system and the bath, i.e. 
Hsb = Yl,j A? ® Bj with j4j G su{7is) and -Bj G su(H B )- In order to protect the 
evolution of S 1 against the effect of the interaction Hsb we seek a cyclic pulse sequence 
as a suitable decoupling interaction. 
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By applying the pulse sequence of an annihilator we get 

H = r^Y J ( U l® 1 ) H ( U 9® 1 ) = 1 S®H B . (14) 
1 1 96G 

This shows that decoupling can be achieved using an annihilator procedure on only 
one of the subsystems (cf . |3l| ) . 

3.4 Inversion of Hamiltonians 

We consider the problem to invert an arbitrary, possibly unknown Hamiltonian, i. e. to 
simulate —H given the Hamiltonian H. 

We can use the following trick |fl^| : by averaging over all elements of G but the 
identity we can invert the Hamiltonian 

U\HU g = -R (15) 

96G\{1} 

because of Lemma || (note that H is traceless). The resulting time overhead is \G\ — 1 = 
d 2 — 1 and the complexity is o? 2 — 1. This can be seen as a generalization of the refocussing 
technique used in NMR. In general, the inverted time evolution will be slower than the 
original one: 

Lemma 2 (Lower bound on inverting) Let r be the greatest eigenvalue and let q 
be the smallest eigenvalue of H. Then r > is a lower bound on the overhead for 
simulating —H by H . 

Proof: Write —H as a positive linear combination of conjugates of H as in Definition |I[ 
Let A m j n (^4) be the smallest eigenvalue of an operator A. Then we have 

i 

The inequality is due to A m i n (j4 + B) > X m i n (A) + A m i n (l?) (see Q, Theorem III. 2) 
for the sum of two Hermitian matrices A and B. Since q is negative it follows that 
—rjq < t. □ 
For the Hamiltonian H = d\&g{d — 1,-1,... ,—1) the overhead is at least d — 1. 
Therefore a lower bound on time overhead for inverting an unknown Hamiltonian is 
d-1. 



4 Universal transformation of Hamiltonians 

In Section ^ we have given a necessary and sufficient condition on the minimal set of 
available control operations in order to enable inversion and cancelling of Hamiltonians. 
If we want to simulate an arbitrary Hamiltonian by any other this condition is not 
sufficient. This can be seen by the following example. Assume that the only control 
operations on C 2 are given by the Pauli-matrices (in their role as unitary operators). If 
the Hamiltonian H := a z is given, conjugation of H by a Pauli-matrix can only lead to 
either H or —H. All the Hamiltonians which can be obtained as average Hamiltonians 
are scalar multiples of H. Hence one cannot simulate e.g. a x . The following concept 
will be useful in order to find groups which enable universal simulation. 
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Definition 4 (Transformer) A subgroup T of SU(d) is called a universal trans- 
former of Hamiltonians iff every M>-linear map L on su(d) (i. e., the set of self-adjoint 
traceless operators) can be written as 

L(A)=Y,PjU}AUj 

j 

with positive real numbers pj and Uj G T. 

The physical of this is that a transformer allows to simulate the Hamiltonian L(H) 
if the unknown Hamiltonian H is present. In |Q] it has been shown that SU(d) is a 
transformer for every dimension d. 

Observation 1 In particular, a transformer is able to simulate an arbitrary Hamilto- 
nian H G su(d) by an arbitrary Hamiltonian H G su(d). 

Remarkably, the condition for a finite group to be a transformer can be characterized 
in terms of irreducibility of certain representations. In contrast to the condition for 
an annihilator, it refers to the adjoint action on the set of operators instead of the 
underlying Hilbert space. 

Definition 5 (Adjoint action) Let G be a finite group and <p a unitary representa- 
tion of degree d, i.e., ip operates onV = C d . We define a linear representation ip^ on 
V ®V by ¥?ad(<?) := fid) ® y>{g) for all g G G, where U denotes complex conjugation of 
a matrix U . We call <p & <\ the adjoint action of ip. Note that this action can be identified 
with the action of G on matrices via conjugation g \— > (M i— > ip(gy M(p(g)). 

In the following we make use of the fact that the algebra generated by the images 
of an irreducible m-dimensional representation i? of a finite group G over the complex 
numbers is equal to the full matrix algebra C mxm . We cite the corresponding theorem 
from ]l3| , Theorem 9.2]. Recall that a representation $ defined over a field F is called 
absolutely irreducible if it remains irreducible when considered over an extension field 
E/F. 

Theorem 2 Let $ be an absolutely irreducible representation of a finite group G which 
has degree m and is defined over the field F . Then 

gdG 

In particular for any m-dimensional irreducible representation over the field C of com- 
plex numbers the vector space generated by the images equals C mxm . 

We now have the necessary prerequisites to characterize finite transformers. 

Theorem 3 (Characterization of finite transformers) 

A finite group T < SU (d) is a transformer if and only if the adjoint representation i? 
given by 

0(i7) := (A ^ U ] AU) 

with U G T acts irreducibly on sl(d) = su(d) + i su{d), i. e., the space of traceless 
operators. 
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Proof: (■<=) Let L be a given R-linear map on su(d) and assume that the adjoint 
action of T is irreducible on sl(d) and denote this representation by $. From Theorem 
|2] follows that the complex linear span of the images of i? is the full matrix algebra acting 
on sl(d). Hence, the mapping L can be written as a complex linear combination of the 
form L:A^J2i QiUjAUi. We now show that the coefficients can be chosen to be real: 
since Xj\AUi is self-adjoint for all i we have L(A) = L(A)^ = ^2i~qlU} AUi. Therefore 
we can write L in the form L : Am J2iPiUjAUi with coefficients pi = ^(qi + £ R. 
Using the inversion scheme of Section |3.4| we can chose the coefficients pi to be positive 
real numbers. 

(=>) Assume that every R-linear map on su{d) can be implemented in the sense of 
Definition |4] using T. Let M be the complex linear span of the maps $(?/) with U £ T. 
The idea is to show that any F £ sl(d), F ^ can be mapped to any other F £ sl(d) 
by a map T £ M. This in turn shows that the adjoint action is irreducible since there 
cannot be a nontrivial invariant subspace. To construct T proceed as follows. Let 
F = Hi + %H 2 with Hi, Hi £ su(d). Assume w.l.o.g. that Hi ^ 0, otherwise multiply 
F by —%. The set Ai contains maps L\ and L 2 with Lx(-ffi) = Hi, Li{H 2 ) = \H\ and 
L 2 (Hi) = H 2 , L 2 (H 2 ) = fiH 2 with n,\ £ R. Then T := Li/(1 + A) +t£ 2 /(l + M) is 
the desired map. □ 



4.1 Finite transformers 

We derive a necessary and sufficient condition for a finite group to be a transformer 
group in the sense of Definition ||. Theorem || shows that the problem to construct 
a finite transformer group is to find for given dimension d > 1 a finite group G and 
an irreducible (unitary) representation (p of G such that the adjoint action becomes 
irreducible if we split off the trivial representation 1 of G. The trivial representation 
is always contained in since the one-dimensional space corresponding to the linear 
span of the identity matrix remains invariant, i. e., c/? a d = l©7r for some representation 
it of G. Abusing the notation we will write <^ a d — 1 to denote the summand ir in this 
decomposition. 

Once we have found a suitable pair (G, ip) with deg(<^) = d this yields a transformer 
group as in Definition 0. For basic results concerning representation theory of finite 



groups we refer the reader to 



Example 2 We examine the case of a two- dimensional system, i. e., d = 2. Starting 
from the Pauli matrices 

/0 1\ / -i \ (I 

a * '= { 1 ) ' *y '= { i ) ' Uz ■= \ -1 

we first note that the group (i ■ a x ,i ■ o~ y ,i ■ a z ) is isomorphic to the quaternion group 
Qs of order 8. This group has an (outer) automorphism of order 3 which permutes the 
Pauli matrices cyclically. This automorphism is given by the matrix 

R:= ' ' / ' ' 



2V" 11 

Setting Sk ■= io~k for k £ {x,y,z} the automorphism is given by i? _1 s z i? = s y , 
R~ l s y R = s z , and R~ 1 s z R = s x . The group generated by the Sk and R is isomor- 
phic to SL(2,F3) ; i. e., the group of 2 x 2 matrices over the finite field F3 which have 
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determinant 1. Let p be the (natural) representation of the matrix group given by 
(s x , Sy, s z , R) . Then the 24 matrices in the image of p form a faithful irreducible 
representation o/SL(2,F3). Choosing the basis {s x ,s y ,s z } of sl(2) we see that the im- 
ages of (^ad — 1 are given explicitly by s x t— > diag(l,— 1 — 1), s y i— ► diag(— 1, 1, — 1), 
s z i — ^ diag(— 1, — 1, 1), and R maps to the permutation matrix corresponding to the 
3-cycle (1,2,3). It is readily verified that this is an irreducible representation. 

Let G be a finite group having an irreducible representation <p such that the images 
of tp are a transformer in the sense of Definition f|. Then necessarily tp must be non- 
monomialQ for otherwise the set of diagonal matrices would be an invariant subspace 
under the action of tp^ — 1. Note that in fact the group SL(2, F3) is the smallest group 
which is not an .M-group, i.e., SL(2,F3) has representations which are not equivalent 
to monomial ones. Therefore a necessary condition for (p to be a transformer has been 
found. 

There is a necessary and sufficient characterization of transformer groups which can 
be verified from the character table alone. Recall that the character x of a representa- 
tion tp is defined by x{d) := tr(<£>(<?)) and that a character is called irreducible iff the 
corresponding representation is irreducible. 

Theorem 4 Let G be a finite group and x be an irreducible character of G with corre- 
sponding representation ip. Then x corresponds to a universal transformer if and only 
if the following identity holds: 

£lxG0| 4 = 2|G|. 

geG 

Proof: The representation p^ — 1 has character values |x(fiO| 2 — 1 f° r all 5 G G since 
tr(<^ad(5 f )) = x(s)x(s)- Recall that the vector space of class functions on G has a 
normalized scalar product given by 

1 1 sec 

for characters Xli X2 of G. A character x is irreducible iff (x\x) = 1- Computing the 
latter scalar product of the character corresponding to ip^ — 1 we obtain 

T^E(i^)i 2 -i)-(ix(5)i 2 -i) 

1 1 sec 



On the other hand this scalar product is equal to 1 due to the irreducibility of p^ — 1. 
Rearranging terms and clearing denominators yields the claimed statement. □ 
In the following we present a transformer for a three dimensional system. The 
minimal group having a representation tp for which tp^ — 1 is irreducible is the linear 

*A representation is called monomial if all representing matrices have the property to contain precisely 
one non- vanishing entry in each row and each column. 



^Ei^)i 4 -^Ei^)i 2 + i 



\G\ 



g&G 
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group GL(3,F2) of invertible 3x3 matrices over the field F2. This is a simple group 
of order 168. As generators of this group we choose the matrices 





1 


0\ 




< 1 


1 




1 


1 









1 




\ 





1 J 









/ 



where x is an element of order 2 and the order of y is 7. The group GL(3, F2) has a 
three-dimensional irreducible representation <p over the complex numbers which on the 
generators x and y is given by the following assignments: 



<p(x) ■= -7= 



2 

7? 



( 



\ 



cos(ff) 
-C 7 4 cos(£) 
-C 7 5 cos(ff) 



-C 7 3 cos(^) 
-cos(ff) 
C 7 cos(ff) 



-C 7 2 cos(ff) 
C 7 6 cos(ff) 



COS^J 



c 7 



c? 



c. 1 



Here • is an abbreviation for and £ 7 denotes the primitive 7-th root of unity e 2m ^ 7 . 
The character of the representation ip takes the values 



3, -1, 1, 0, C7 + C7 + C 



7- 



c 7 3 + c 7 5 + c 7 6 



on the conjugacy classes of GL(3, F2). Consulting the character table of GL(3,F2) in 
the Atlas || p. 3] we find that ip is irreducible. The representation ip a ^ — l has character 
values 



8, 0, 0, -1, 1, 1 

from which follows that it is also irreducible, again by checking the character table of 
GL(3, F2). Overall we obtain that the representation tp of GL(3, F2) yields a transformer 
of size 168. Using the Neubiiser catalogue used in MAGMA and GAP, cf. [f|, 3C] we 



performed an exhaustive search over all groups of smaller sizes which has shown that 
this indeed is the minimal possible group size. 

In Tabular [l] we summarize the results of this search. Groups of sizes up to 255 have 
been considered. The number in the Neubiiser catalogue is given such that for instance 
the first row of this table corresponds to the group (in GAP syntax) SmallGroup (24 , 3) 
which has been studied in Example ^. Note that we only give transformer groups which 
act faithfully. 



4.2 Lower bound on the overhead 

In the following we derive a lower bound on the time overhead for simulating a Hamilto- 
nian using an arbitrary transformer. We need some results on majorization and doubly 
stochastic matrices (cf. |^]|] for a summary). Let x = {x\, . . . , Xj) and y = (yi, . . . , y^) 
be two d-dimensional real vectors. We introduce the notation j to denote the com- 
ponents of a vector rearranged into non-increasing order, so 2^ = (x[, . . . ,xi), where 
(x| > X2 > • • • > Xj). We say that x is majorized by y and write x -< y, if 

k k 

».;•:>>.;• 

j=i 3=1 
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Group size 


Numbers in library 


Dimension 


24 


3 


2 


48 


28, 29, 33 


2 


72 


3, 25 


2 


96 


67, 74, 192 


2 


120 


5 


2 


144 


36, 121, 122, 157 


2 


168 


22 


2 


168 


42 


3 


192 


187, 204, 963 


2 


216 


3, 38 


2 


216 


88 


3 


240 


93, 102, 103, 154 


2 



Table 1: Transformer groups of small sizes 



for k = 1, . . . , d — 1, and with equality when k = d Q. 

Let Spec(AT) denote the spectrum of the hermitian matrix X, i.e. the vector of 
eigenvalues. Ky Fan's maximum principle gives rise to a useful constraint on the 
eigenvalues of a sum of two Hermitian matrices C := A + B, that 

Spec(^ + B) ~< Spec(A) + Spec(£) . (16) 

This permits us to derive a lower bound on the simulation overhead. 

Lemma 3 (Lower bound) A lower bound on the overhead of simulating H by H is 
given by the minimal r such that 

Spec(H) ~< rSpec(tf) . (17) 

Proof: This follows from Definition |l| and inequality ((l^) . □ 
We now consider the question when this lower bound can be attained. Let H = 
^^ =1 Aj|«)(i| and H = ^f=i/ i «IO(*l where \i) is a basis of eigenvectors. Let r be 
minimal such that /2 -< rA. We set d' := d — 1. We have . . . , Hd>) ~< (tAi, . . . , rA^) 
since Y2i tM = Y2i Ai = (the Hamiltonians are traceless) . This is equivalent to the 
existence of a doubly stochastic matrix D with -D(Ai, . . . , A^-i) T = (fi\, ■ ■ ■ , Hd-i) T ■ 
By Birkhoff 's theorem we can decompose a doubly stochastic matrix as a convex sum 
of permutations, i. e. 

D = Y J PaU a (18) 
o-es 

where U a is the permutation matrix associated to a, i.e. maps the basis vector \i) 
to \a(i)), and and S is a subset of the symmetric group Sd'- The d! x d' doubly 
stochastic matrices form a (d' 2 — 2d' + l)-dimensional convex set. The extreme points 
are the permutation matrices. Caratheodory's theorem guarantees that a point in a 
m-dimensional compact convex set may be expressed as a convex combination of at 
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most m + 1 extremal points of that set. Therefore every doubly stochastic matrix can 
be written as a convex combination of at most d' 2 — 2d' + 2 = (d— 2) 2 + 1 permutations. 
We view the matrices U a as d x d matrices that fix the basis vector | d) . We have 

H = rY J P*UlHU a (19) 

The lower bound can be attained in particular if the transformer contains the matrices 
permuting the eigenvectors of H and and the matrix realizing the base change between 
the eigenvector basis of H and H. 

Let H be the system Hamiltonian. Let a be the cyclic shift, i. e. a(i) = i + 1 mod d 
and U a the corresponding matrix. Than we can switch off the Hamiltonian 

E^>^=0 (20) 

3=0 

with complexity at most d provided that we can perform the shift with the transformer. 
In this case we can also invert it with complexity and overhead at most d—1. However, 
when the Hamiltonian is not known we need at least d 2 operations to switch it off. 



5 Simulation of bipartite Hamiltonians 

The following theorem shows that all bipartite Hamiltonians can be simulated by any 
Hamiltonian (with non-trivial coupling and non-trivial local terms) provided that the 
set of available unitary transformations contains a transformer for each of the subsys- 
tems. Let B = {a a \ a = 1, . . . , d 2 — 1} be a basis of su(d). 

Theorem 5 Let an arbitrary interaction 

H := ^2 Jap&a ® crp + a®l + l®b 

af3 

be given with a,b £ su(d) with a, b ^ 0. Let T\ and T2 be transformers of the left and 
the right subsystem, respectively. Assume that it is possible to implement all unitary 
transformations of the form U <8> V with U G T\ and V G T2. Then H can be used for 
simulating any arbitrary H with H := J2 a /3 Ja/3°~a ^c^ + a^l-l-l^fe i. e., there are 
positive numbers Tj, Uj £ t\ and Vj £ T2 such that 

k = Y. t M ® V 1) H ( U 3 ® V,) ■ (21) 

3 

Proof: We first consider the case that the local terms are all zero. Write H in the 
form 

H = Y,A j ®B j , 

3 

where Aj and Bj are elements of su(d) and all Bj are linearly independent and all 
Aj are nonzero. Then H can be transformed into any interaction of the form C D 
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with arbitrary C, D £ su(d). This can be done by choosing M-linear maps L\ and 
L 2 on su(d) with L\{A\) = C, L 2 {B X ) = D, and L 2 (Bj) = for j ^ 1. Since 7i 
and 72 are universal transformers one can find positive numbers dj and fj and unitary 
transformations Uj S 71 and Vj £ T2 such that 

L!(.) = ^ • and L 2 (-) = £ ./'l'; • 1} . 

j j 

Hence we obtain 

X] ® ^ ® ^) = c ® ^ ■ 

This proves that we can simulate each tensor product operator C ® D. By setting 



C '■= Jap&a an d D := op we can simulate each term in eq. (21). Hence it is possible 
to simulate H. Let H contain local terms. Starting from H we can simulate the 
Hamiltonian 

H' = J a pa a 8) ap + a' (g) 1 + 1 <g) 6' , 

a/3 

that coincides with the desired Hamiltonian 77 except for the local terms. Starting 
from H we can simulate the Hamiltonian (a — a') ® 1 by applying an annihilator on 
the right and suitable transformations on the left. Finally, we use a similar scheme for 
l(g> (&-&')• D 



6 Conclusions 

We have shown that there are finite groups of unitary control operations which allow 
to simulate an arbitrary Hamiltonian H by another arbitrary H, i.e., a system with 
Hamiltonian H can be driven to evolve as if its Hamiltonian were H. This can be 
accomplished using fast sequences of control operations interspersing the natural time 
evolution. We even found finite groups which allow to solve the following more general 
control problem: for every linear trace preserving map L on the set of self-adjoint op- 
erators the system can be made to evolve according to the Hamiltonian L(H) although 
its true unknown Hamiltonian is H. We have called such groups transformer groups 
and showed that a finite group G has this property if and only if its adjoint action 
on the set of traceless operators is irreducible. We have characterized finite groups 
G with this property using characters of representations of G on the Hilbert space. 
This criterion allows to perform an exhaustive search over groups of small order (up to 
255). We found transformer groups for two and three dimensional quantum systems. 
It remains an open problem to construct finite transformer groups for all dimensions. 
In bipartite systems, every non-trivial interaction can simulate any other provided that 
transformer groups on each subsystem can be implemented. 
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